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A new method of accessing information on the symmetry free energy from yields of fragments 
produced in Fermi-energy heavy-ion collisions is proposed. Furthermore, by means of quantum 
fluctuation analysis techniques, correlations between extracted symmetry free-energy coefficients 
with temperature and density were studied. The obtained results are consistent with those of 
commonly used isoscaling techniques. 
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I. INTRODUCTION 

The fundamental goal of studying strongly excited nu¬ 
clei is to characterize the equation of state (EOS) of nu¬ 
clear matter over a wide range of temperature, density, 
pressure and isospin. This quest to characterize nuclear 
phenomena at intermediate energies is still a matter of 
intensive investigations. The characterization of the EOS 
of nuclear matter plays a key role for various phenomena 
in nuclear astrophysics, nuclear structure, and nuclear 
reactions [1-5]. The EOS for asymmetric nuclear mat¬ 
ter is usually expressed as a term related to symmetric 
matter and a term which takes into account the isospin 
asymmetry of the system. The latter is referred to as 
the symmetry energy. In early studies, the symmetry en¬ 
ergy coefficient (E sym ) of nuclei was extracted by fitting 
the binding energy in their ground state with various 
versions of the liquid drop mass formula. The proper¬ 
ties of nuclear matter are afterwards determined by the¬ 
oretically extrapolating the nuclear models designed to 
study the structure of real nuclei which are cold (T = 0), 
nearly symmetric (N ~ Z) and at the saturation density 
(/?o ps 0.16 fm -3 ). In contrast to the value of E sym at p 0 
and T = 0, the behavior of E sym with temperature (T) 
and density (p) is still being mapped out. 

Many experimental and theoretical investigations have 
been devoted, in recent years, to estimate the behavior 
of E sym as a function of T and p. Among these efforts, 
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measurements of the giant dipole [6], pygmy dipole [7] 
and giant monopolc [8] resonances in neutron-rich nu¬ 
clei, neutron-proton emission ratios [9], isospin diffusion 
[10], collective flows [11] and fragment isotopic ratios [12- 
15] have provided constraints on the density dependence 
of the symmetry energy at subsaturation densities. Re¬ 
cently, a large body of experimental data from studies 
of heavy-ion collisions [16-19] have been used to extract 
the free symmetry energy and the symmetry energy at 
subsaturation densities and moderate temperatures. In 
those studies, isoscaling parameters deduced from iso¬ 
topic yields measured in two similar reactions with dif¬ 
ferent isotopic composition were used to access the sym¬ 
metry free energy coefficients. The symmetry energy 
coefficients were in turn derived using model calculated 
symmetry entropies together with experimental symme¬ 
try free energy coefficients. 

To a good approximation, the EOS of asymmetric nu¬ 
clear matter can be written as 

E(p, T, m) = E(p , T,m = 0) + E sym (p , T)m 2 + 0(m 4 ) , 

(1) 

where p = p p + p n , m = (N — Z)/A is the neutron-proton 
asymmetry and E(p, T,m = 0) is the EOS of symmetric 
nuclear matter [20, 21]. As a general representation of 
the symmetry energy coefficient, the following definition 
has been considered 

E sym (p, T) = 1 [E(p, T, 1) + E(p, T, -1)] - E(p, T, 0) . 

( 2 ) 

We recall here that in most cases, the symmetry energy 
is connected to the isotopic yields through the relation¬ 
ship 



where C sym is the symmetry energy coefficient, Z\, A\ 
and Zi, A 2 are respectively the charge and the mass num- 
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bers of system 1 and 2 (system 2 being richer in neutrons 
than system 1), and T is the common temperature of 
the two systems [13, 15, 22]. However, at sufficiently low 
densities, C sym was shown to be substituted for the sym¬ 
metry free energy ( F sym ) and is related to E sym through 
E sym = F sym + TS S ym, with S sym the symmetry entropy 
[16-19]. This implies that at low densities the symmetry 
entropy contribution to E sym becomes significant as clus¬ 
tering increases the binding energy and therefore reduces 
the entropy in symmetric matter [23]. 


fragmenting source are determined using the quantum- 
fluctuation method, fully described in Refs. [37-41], This 
is the first time that experimental fragment yield data 
analyzed within the Landau free-energy framework are 
used to determine symmetry free energy coefficients in a 
completely self-consistent manner. 


III. EXPERIMENTAL DETAILS AND EVENT 
SELECTION 


II. FORMALISM 

Heavy-ion collisions in the Fermi-energy domain are 
dominated by nuclear fragmentation and their studies 
provide information about the properties of nuclear mat¬ 
ter at moderate temperatures and sub-saturation den¬ 
sities. Several studies [24-26] have shown this energy 
domain to be the region of the nuclear liquid-gas phase 
transition. In our recent works [27-32] , we have analyzed 
fragment yield data to investigate the nuclear phase tran¬ 
sition using the Landau free-energy approach [33, 34]. 
In such an approach, the key assumption is that in the 
vicinity of the critical point, the fragment free energy per 
nucleon ( F ) relative to the system temperature (T) can 
be expanded in a power series in the fragment’s neutron- 
proton asymmetry m as given by the relation 


F 1 2 1, 4 1 6 H 

— = -am + -bm ± -cm b - — m , (4) 

where m = (Nf — Zf)/Af , and Nf, Zf, and Af are the 
neutron, proton, and mass numbers of the fragment, re¬ 
spectively. The quantity m behaves as an order param¬ 
eter, H is its conjugate variable and the coefficients a, 
b, and c are fitting parameters. We recall that based 
on a modified Fisher model [27, 35, 36], fragment yields 
are proportional to Aj T e~^ F ' T ' )A f near the critical point; 
with r as the critical exponent. 

From Landau’s free energy equation, using an anal¬ 
ogous expression of Eq. 2 for F sym (composed of pure 
neutron, pure proton, and symmetric nuclear matter), 
one can obtain the following expression: 


TP 

x sym 

T 



( 5 ) 


The parameters of the Landau’s equation are related to 
the state variables of the fragmenting system and have 
been shown to depend on its proton-neutron concentra¬ 
tion and excitation energy [30, 32]. This suggests that 
these parameters could be used to directly obtain infor¬ 
mation about the symmetry free energy which is a com¬ 
ponent of the nuclear EOS. 

In this paper, we report on experimental symmetry 
free energy coefficients extracted using the Landau free- 
energy approach. The temperature and density of the 


The experiment was performed at the K-500 super¬ 
conducting cyclotron facility at Texas A&M University. 
Beams of 64 Zn, 70 Zn and 64 Ni at 35 MeV/A were used 
to respectively irradiate 64 Zn, 70 Zn and 64 Ni targets. 
The 47 t NIMROD-ISiS array [42, 43] was used to col¬ 
lect charged particles and free neutrons produced in the 
reactions. More details of the experiment have been 
given in Refs. [44-46]. An excellent energy resolution was 
achieved, allowing isotopic resolution of charged particles 
up to Z = 17 and elemental resolution up to the charge 
of the beam. 

We reconstruct the primary hot nuclear system in or¬ 
der to select events of similar character. For this anal¬ 
ysis, we are interested in studying equilibrated systems, 
so events are selected in the following way: Fragments 
that do not originate from an equilibrated quasiprojec¬ 
tile (QP) source were excluded with the condition that 
the longitudinal velocity of fragments with Z = 1, 2, > 3 
be in the range of ±65%, ±60%, and ±40%, respectively, 
of the velocity of the heaviest fragment in the event. The 
sum Z is selected from 21 to 30, the sum A of fragments is 
selected to be 54 to 64. To select roughly spherical events, 
which are spatially equilibrated, the quadrupole moment 
Qshape was required to be -0.3 < log 10 (Q shape) < 0.3 
where Qshape = J2 pI/J 2 \Pt > an d Pz and p t are respec¬ 
tively the longitudinal and transverse momenta of the 
fragments comprising the QP. 

The QP system was reconstructed from events in which 
all charged particles and free neutrons were isotopically 
identified. The neutron ball provided event-by-event ex¬ 
perimental information on the free neutrons emitted dur¬ 
ing a reaction. The number of free neutrons emitted by 
the QP was deduced from the total measured number of 
neutrons, background, and efficiencies for measuring neu¬ 
trons produced from QP and quasitarget sources [47, 48] . 
In order to minimize contributions from collective effects, 
which are predominant in the beam direction, the exci¬ 
tation energy of the reconstructed QPs was estimated 
from transverse kinetic energy of the charged particles, 
the neutron multiplicity, the average neutron kinetic en¬ 
ergy determined using the Coulomb-shifted proton en¬ 
ergy distribution, and the energy needed for the breakup 
(Q- value). A detailed description of this method of re¬ 
construction is given in Refs. [47, 49]. Events were sorted 
in 8 excitation energy bins, 1 MeV/A wide, from 2.5 to 
9.5 MeV/A 
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IV. TEMPERATURE AND DENSITY 

The temperatures of reconstructed QP sources are 
obtained with the quadrupole momentum fluctuation 
method. The method is described in a very detailed 
way in [37-41, 50]; here we only briefly outline it. The 
quadrupole momentum is defined as Q xy = p x — p y us¬ 
ing the transverse components p x and p y of the particle’s 
momentum in the frame of the QP source. If the correct 
quantum distribution for fermions is used, the variance 
of Q xy is related to the temperature by 

(vly) = 4m lartT 2 F Q c , (6) 

where m par t is the mass of the particle being used as the 
probe and Fqq is the quantum-correction factor. The 
quadrupole is defined in the transverse direction in order 
to minimize nonequilibrium effects which may manifest 
in the beam direction. Equation 6 was solved numerically 
for a Fermi gas and the quantity Fqq was parametrized 
in terms of the temperature relative to the Fermi energy 
( [T/ef ) which in turn was parametrized in terms of the 
normalized multiplicity fluctuation. The nucleon den¬ 
sity p is therefore determined from the Fermi-energy re¬ 
lation Ef = ef 0 (p/po) 2 ^ 3 with £f 0 and pa respectively the 
ground-state values of Fermi energy and nucleon density. 
Coulomb corrections were later applied to derived tem¬ 
peratures and densities. This is done by a method bor¬ 
rowed from electron scattering where the Coulomb field 
is taken to be the Fourier transform of the Coulomb po¬ 
tential of the source. In such an approach, as described 
in Refs. [40, 41], the equations of quadrupole momen¬ 
tum fluctuation, the average multiplicity, as well as the 
multiplicity fluctuation which contain the Coulomb field 
term are numerically solved to derive the temperature 
T, the density p, and the volume V of the system. This 
method can be considered reliable as from model calcula¬ 
tions, temperatures of protons after Coulomb corrections 
were similar to those of neutrons. In Ref. [51], where we 
have applied this method to experimental data, it was 
observed that these Coulomb corrections lower temper¬ 
ature values by almost 2 MeV while their effects were 
small on derived densities. The error in applying the 
Coulomb corrections arises from the uncertainty in the 
source charge. We varied the source charge by ± 2 units 
and the estimated errors are respectively ±2% for the 
densities and ±6% for the temperatures. 

Previously [31, 32, 50] the data were sorted into four 
different QP asymmetry bins of width 0.05, ranging from 
0.04 to 0.24. In a subsequent paper [51], it was shown 
that values of T and p that correspond to an asymme¬ 
try bin width close to zero (as it should be for fixed A 
and Z) could be obtained by averaging values for all four 
asymmetry bins. Therefore, in this paper we have re¬ 
ported averaged values of temperatures and densities for 
all four asymmetry bins. On the same grounds, F sym 
values presented here are also obtained from averaged 
values of Landau’s equation fitting parameters from the 
four asymmetry bins. 


V. RESULTS AND DISCUSSION 

Figure 1 shows the free energy ( F/T ) values as a func¬ 
tion of fragment’s neutron-proton asymmetry m at an ex¬ 
citation energy of 5.5 MeV/A of the QP. As the efficiency 
for measuring neutrons differs from the efficiency for mea¬ 
suring charged particles, only charged-particle yields are 
used in the Landau’s equation fitting. The form of the 
fitting function has physical restrictions on it. Since the 
symmetry energy is isospin symmetric, only even powers 
of m appear, aside from the external field. Furthermore, 
at extreme values of asymmetry, the free energy must 
not be decreasing toward negative infinity. Within these 
constraints, the form must also have a minimum at zero 
as dictated by the data. A parabolic fit describes the 
data only around zero, but misses entirely the points at 
both to = — 1 and to = 0.5. Our ability to measure 
this point is very useful to constrain F/T. Having ex¬ 
cluded the quadratic fit, we next rule out all quartic fits 
immediately on the grounds that the function must be 
rising at extreme values of m and still have a minimum 
at 0. A sixth order (in even terms only) is the next sim¬ 
plest polynomial that satisfies the physical constraints. 
The values of the free energy obtained were corrected 
for pairing [27, 32], similar to the mass formula, and a 
good scaling is seen in the figure. The solid line (Landau 
Fitl) is a fit to the data with all parameters as free pa¬ 
rameters. In Ref. [32], the parameter c was observed to 
be almost constant, within uncertainties, over the entire 
range of the QP excitation energy. Here we additionally 
fitted the data fixing c = 115. This is represented by the 
dashed line (Landau Fit2). It is observed that the two 
fitting curves provide a good fit to the free-energy data. 
However, at the two extreme minima (at large to values) 
where we have no data points the two curves are slightly 
different. The values of a, b and c corresponding to the 
solid line were obtained as 15.527±0.041, -91.786±0.484 
and 100.81±0.615, respectively. After fixing c = 115 
(dashed line), the values of a and b were obtained as 
16.289±0.024 and -102.871±0.058, respectively. This re¬ 
sulted in similar values of F sym /T, calculated using Eq. 5, 
as 1.619±0.16 and 1.593±0.019. The values of parame¬ 
ters a and b of Landau’s equation used in the rest of the 
discussion were obtained by fixing c = 115. In this way, 
estimated errors on extracted F sym values were signifi¬ 
cantly minimized. The appearance of the three minima 
is a signature of a first-order phase transition of the sys¬ 
tem [27-32]. 

As there are no experimental data points in the re¬ 
gion of the two minima of the fit and one must rely on 
the proton point at to = — 1, we have investigated the 
uncertainty on the F sym values introduced by the fit pa¬ 
rameters if this point was in error by some amount. A 
15% variation of the F/T value at to = —1 resulted in 
a 19% change of the F syrn value. We estimated the sys¬ 
tematic uncertainty caused by the isolated proton point 
on the F sym values to be less than 20%. 

The symmetry free energy coefficients determined from 
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FIG. 1: (Color online) F/T values as a function of fragment’s 
neutron-proton asymmetry m for an excitation energy of 5.5 
MeV/A of the QP. The solid line (Landau Fitl) is a fit to the 
data with all Landau’s equation parameters as free parame¬ 
ters while the dashed line (Landau Fit2) represents a fit to 
data fixing c = 115. Error bars corresponding to statistical 
errors are smaller than the symbols. 



FIG. 2: (Color online) (a) Symmetry free energy vs temper¬ 
ature. (b) Symmetry free energy vs density. The symmetry 
free energy coefficients ( F syrn ) are extracted from Landau’s 
free energy approach. Temperatures and densities are de¬ 
rived from the quantum-fluctuation method with protons as 
the probe particle. Full circles and squares correspond respec¬ 
tively to results without and with Coulomb correction. For 
comparison, F syrn values obtained from isoscaling techniques 
along with T and p from Albergo [17] and coalescence [18] 
methods are also plotted. Statistical errors are indicated by 
the bars and are not shown when smaller than the symbols. 


the Landau free-energy technique are displayed as a func¬ 
tion of the system temperature T (Fig. 2(a)) and den¬ 
sity p (Fig. 2(b)). T and p for this work were derived 
from the quantum-fluctuation method without and with 
Coulomb corrections, while Fsyrn values are derived from 
the Landau equation parameters. The present data set is 
compared to previously published work, which used the 
Albergo method. The Albergo method is a double ratio 
technique that evaluates the temperature and density of 
equilibrated nuclear regions using the yields of different 
light nuclides (d, t, h, a). Application of this technique 
assumes that thermal equilibration and chemical equili¬ 
bration have been attained [16, 52]. Data represented by 
inverted triangle symbols (Albergo) are taken directly 
from Ref. [17], where T and p were determined. The 
isoscaling method was used to extract F sym values. A 
comparison is also made to previous work that extracted 
T and p using a combination of the Albergo method and a 
coalescence model. This work is shown in full up-triangle 
symbols, the points are taken directly from Ref. [18]. The 
overlap of the errors from this work reflect that the bin¬ 
ning of the data was not done as a function of tempera¬ 
ture, density, or free symmetry energy, but of fragment 
surface velocity, which is correlated to the emission time. 
The effect of this is seen in the size and proximity of 
the error bars. In the coalescence model, the momen¬ 
tum space densities of ejected light composite particles 
are directly related to those of the ejected nucleons of 
the same velocity. The phase space correlations, which 
lead to cluster formation, may therefore be parametrized 
in terms of the radius of the momentum space volume, 
Pq, within which the correlations exist. In this way the 
double isotope yield ratio at equal velocity is used to de¬ 
termine the temperature. From the relationship between 
the coalescence parameter Pq and the volume of the emit¬ 
ting system, the density is derived [53, 54] . Isoscaling was 
also used to extract F sym in this case. 

It is observed that there is a fair agreement be¬ 
tween temperature-dependent F sym results. However, 
the density-dependent F syrn results with Coulomb cor¬ 
rections significantly deviate from those obtained without 
Coulomb corrections and Coalescence and Albergo data 
as well. Neverthless, it is amazing to see that within 
the error bars temperature- and density-dependent sym¬ 
metry free energies derived from different methods agree 
with each other to a remarkable degree. In Ref. [18], 
where the colescence data have been taken, the quoted er¬ 
rors on the temperatures are 10% at low density evolving 
to 15% at the higher densities. The error in the deriva¬ 
tion of the density was estimated to be in the order of 
17%. 

From the values of F symi T and p , we examine in Fig. 3 
the symmetry free-energy density (3 sym = F sym x p) 
against T. It is observed that Zsym monotonically in¬ 
creases as T increases and some of the differences seen 
between curves in both panels of Fig. 2 are less evident 
except for the curve represented by full squares where 
Coulomb corrections have been accounted for. There- 
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FIG. 3: (Color online) Symmetry free-energy density (3sym = 
Fsym x p) as a function of the system temperature T. Er¬ 
ror bars represent statistical errors and are not shown when 
smaller than the symbols. 


fore, plots of F sym values as a function of T{p) could 
be misleading, since T and p vary simultaneously. The 
Zsym as a function of temperature makes use of all the 
three quantities that are experimentally accessed. This 
quantity displays a clear deviation among the different 
methods especially when dealing with the Coulomb cor¬ 
rection, suggesting some model dependence in dealing 
with the Coulomb term. This ambiguity could be fur¬ 
ther studied by deriving density and temperature from 
neutron fluctuations as well. We expect from model cal¬ 
culations [41] that quantities derived using the protons 
should be similar to those of neutrons after the Coulomb 
correction. If this is not true, then further work is needed 
to understand the role of Coulomb. In the coalescence 
model a Coulomb correction is applied through a shift of 


the measured kinetic energy spectra of the particles of 
interest. From such a Coulomb shift, knowing the charge 
and mass of the emitting system, it is possible to derive 
its density [55] which can be compared to the obtained 
values from coalescence [18, 19]. If the densities obtained 
do not agree with each other, then a simple shift of the 
energy distribution is not sufficient to describe the role 
of Coulomb. 


VI. CONCLUSIONS 


In conclusion, we have shown that symmetry free en¬ 
ergy coefficients of nuclear systems can be extracted from 
fragment yield data produced in Fermi-energy heavy-ion 
collisions by employing Landau’s free-energy approach. 
The temperature- and density-dependent symmetry free 
energies have been observed to be consistent with those 
derived from isoscaling analyses. This is the first time ex¬ 
perimental fragment yield data analyzed within the Lan¬ 
dau description have been used to determine symmetry 
free energy coefficients. We have found some discrepancy 
among different methods possibly because of the different 
handling of Coulomb corrections. Precise measurement 
of the neutron distribution function might help to solve 
this ambiguity. The estimation of entropic contributions 
to the symmetry free energy in order to derive symmetry 
energy coefficients is currently being given special atten¬ 
tion. 
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